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Abstract 
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Q ■ The finite temperature Brueckner-Hartree-Fock approach is extended by in- 

' troducing a microscopic three-body force. In the framework of the extended 

(N ■ 

model, the equation of state of hot asymmetric nuclear matter and its isospin 
\ dependence have been investigated. The critical temperature of liquid-gas 

' phase transition for symmetric nuclear matter has been calculated and com- 

CN \ pared with other predictions. It turns out that the three-body force gives 

■ a repulsive contribution to the equation of state which is stronger at higher 

density and as a consequence reduces the critical temperature of liquid-gas 
phase transition. The calculated energy per nucleon of hot asymmetric nu- 
^ , clear matter is shown to satisfy a simple quadratic dependence on asymmetric 



parameter j3 as in the zero-temperature case. The symmetry energy and its 



^ , density dependence have been obtained and discussed. Our results show that 

' the three-body force affects strongly the high-density behavior of the symme- 

try energy and makes the symmetry energy more sensitive to the variation 
of temperature. The temperature dependence and the isospin dependence 
of other physical quantities, such as the proton and neutron single particle 
potentials and effective masses are also studied. Due to the additional repul- 
sion produced by the three-body force contribution, the proton and neutron 
single particle potentials are correspondingly enhanced as similar to the zero- 
temperature case. 
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I. INTRODUCTION 



One of the most important motivations of heavy ion physics is to investigate the prop- 
erties of nuclear matter under extreme conditions of temperature, density and isospin. The 
equation of state (EOS) of asymmetric nuclear matter at finite temperature plays an im- 
portant role in understanding the dynamics of heavy ion collisions at intermediate and high 
energies. Theoretical simulations [1] indicate that the main reaction dynamics and observ- 
ables including collective flow, balance energy, isospin equilibrium, pre-equilibrium nucleon 
emission, isotopic scaling and isospin diffusion are quite sensitive to the isospin dependent 
part of the nuclear EOS, especially to the density dependence of the symmetry energy. 
The properties of hot and dense asymmetric nuclear matter are also of great interest in 
astrophysics especially in connection with the dynamics of supernova explosions and the 
thermal evolution of neutron stars [2] . The temperature and isospin dependence of the EOS 
of asymmetric nuclear matter is crucial for a reliable description of the thermodynamical 
evolution and the structure of the newborn neutron star ( protoneutron star ) formed in 
the latest stage of a type-11 supernova collapse [2]. Motivated by the general significance in 
heavy ion physics and astrophysics, the thermodynamical properties of nuclear matter has 
been investigated extensively within various theoretical models such as the non-relativistic 
Hartree-Fock approach [3] and the relativistic mean field theory (RMT) [4,5]. 

Due to the Van der Waals nature of the nucleon-nucleon (NN) interaction, it is generally 
expected that nuclear matter is likely to exhibit a liquid-gas phase transition [6]. Theoreti- 
cally, considerable effort has been devoted to establishing the EOS of hot symmetric nuclear 
matter and discussing the critical phenomena [7-13]. Although almost all the predicted 
EOSs of infinite symmetric nuclear matter display a typical Van der Waals behavior, the 
obtained values of the critical temperature Tq for the phase transition are distributed in 
a wide range from 8MeV to 20MeV. Within the non-relativistic framework, the calculated 
value of Tc for infinite symmetric nuclear matter is 15 — 20MeV [9,10] depending on the 
choice of the theoretical models and the NN interactions. The RMT theory gives a critical 
temperature of about 14MeV [5] . The value of Tc predicted from the Dirac-Brueckner ( DB 
) approach is as low as ~ 8MeV in Ref. [12] and about lOMeV in Ref. [13]. 

On a microscopic basis, the properties of hot asymmetric nuclear matter have been 
studied in Ref. [14] in the framework of the finite temperature Brucckncr-Hartrcc-Fock 
( FTBHF ) approach using the separable version of the Paris two-body force. As well 
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known, three-body forces are necessary for reproducing the empirical saturation properties 
of cold symmetric nuclear matter in a non-relativistic microscopic approach [15-19]. Two 
kinds of three-body force ( TBF ) have been adopted in the Brueckner-Hartree-Fock ( BHF 
) formalism [15-18]. One is the semi-phenomenological TBF [20] which has two or few 
adjustable parameters determined by fitting the empirical saturation density and energy of 
the cold symmetric nuclear matter in the BHF calculations [18]. The other is the microscopic 
TBF based on meson exchange coupled to the intermediate virtual excitations of nucleon- 
antinucleon pairs and nucleon resonances [15-17]. 

In the present paper, our aim is to extend the FTBHF approach by introducing the 
microscopic TBF and to investigate the properties of hot asymmetric nuclear matter. Spe- 
cial attention has been payed on the effects of the TBF contribution and the temperature 
dependence. The two-body reahstic NN interaction adopted in the present calculations is 
the Argonne Vis ( ^Vls ) potential [21]. Our calculations indicate that the TBF affects 
not only the EOS of hot asymmetric nuclear matter and the critical temperature of the 
liquid-gas phase transition, but also the single-particle ( s.p. ) properties such as the proton 
and neutron s.p. potentials and effective masses. The present paper is organized as follows. 
We will describe briefly the model in Section II. The numerical results are presented and 
discussed in Section III. In Section IV a summary of the present work is given. 



II. THEORETICAL MODEL 



In general, three parameters are required to specify a given thermodynamical state of hot 
asymmetric nuclear matter, i.e., the total nucleon number density p, the isospin asymmetry 
parameter /3 — BjizEe ^nd the temperature T. At zero-temperature, the proton and the 
neutron Fermi momenta are given by kp — [Stt^p,-]^/^ with r — p or t — n, where p„ = | 
{1 + P)p and pp — ^{1 — P)p are the proton and neutron number densities, respectively. The 
formalism of the Brueckner-Bethe-Goldstone (BBG) theory for cold asymmetric nuclear 
matter can be found in Ref. [22]. The extension to finite temperature is given in Ref. 
[14]. In the following, we give a brief review for completeness. The starting point of the 
BBG approach is the Brueckner reaction matrix G which satisfies the following generalized 
Bethe-Goldstone equation, 

GrAP, /3, T, c) = ^ + . E ^^'^i?7'^^'^^,\ GrAP, T, CO) (1) 
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where v — V2 + V^^ is the nucleon-nucleon (NN) interaction and uj is the starting energy. 
In the present calculations the Argonne Vig (AVis) potential [21] is adopted as the bare 
two-body force V2 and V^^*'^ is the TBF contribution given by Eq.(7)(see below). The finite 
temperature Pauli operator Qr,T' is simply an extension of the zero-temperature one, i.e, 



Qr,r' = QrAkl,k2,P,P,T) = [1 - fr{h, p, (3,T)][l - U^k^, p, (5 ,T)\ 

At finite temperature the Fermi distribution is expressed as, 

' er{k) - p,r 



(2) 



fr(k,p,P,T)^ 



1 -I- exp 



T 



-1 



(3) 



where fj,r — Pr{k,p,T) refers to the chemical potential and er{k) is the s.p. energy. For 
any given density and temperature, we can calculate the chemical potential Pr from the 
following implicit equation self-consistently by iteration. 



Pr-YrY.fr{k,P.P.T) = W 
^ k ^ k 



1 -|- exp 



'er{k) - Pr 



(4) 



The s.p. energy is given by. 



er{k) = er{k,p,/3,T) 



2^,2 



h'k 

2m 



+ Ur{k,p,P,T) 



(5) 



where Ur{k, p, (3,T) is the proton or neutron s.p. potential. Within the FTBHF framework 
the s.p. potential is calculated from the on-shell anti-symmetrized G matrix, 

Ur{k,p,(3,T) = ^J2T.fr'ik',p,P,T){kk' I Gr,r'{p,f3,T,er{k) + er^{k')) I kk')A (6) 

■r' k' 

In the present calculations, we adopt the continuous choice [23] for the s.p. potential. On 
the one hand, the continuous choice has been shown to provide a much faster convergence 
of the hole-line expansion for the energy per nucleon in nuclear matter at T = than the 
gap choice [24]. On the other hand, in the continuous choice, the s.p. potential describes 
physically the nuclear mean field felt by a nucleon in nuclear medium. 

The microscopic TBF adopted in the present calculations is constructed from the meson- 
exchange current approach [15]. Its components are depicted diagrammatically in Fig. 1, 
taken from Ref. [15]. Four important mesons tt, p, a and lo are considered [19]. The TBF 
contains the contribution of the two-meson exchange part of the NN interaction medium- 
modified by the intermediate virtual excitation of nucleon resonances, the term associated 
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to the non-linear meson-nucleon coupling required by the chiral symmetry, the simplest con- 
tribution rising from meson-meson interaction and finally, the two-meson exchange diagram 
with the virtual excitations of nucleon-antinucleon pairs. The meson masses in the TBF have 
been fixed at their physical values except for the virtual cr-meson mass which has been fixed 
at 540MeV according to Ref. [15]. This value has been checked to satisfactorily reproduce 
the AVis interaction from the one-boson-exchange potential (OBEP) model [17]. The other 
parameters of the TBF, i.e., the coupling constants and the regular masses (form factors), 
have been determined from the OBEP model to meet the self-consistent requirement with 
the adopted AVis two-body force. For a more detailed description of the model and of the 
approximations we refer to Refs. [15-17]. 

The TBF effect is included in the self-consistent Brueckner procedure along a similar way 
to the zero temperature case [15,16], where an effective two-body interaction is constructed to 
avoid the full three-body problem. A detailed description and justification of the procedure 
can be found in Ref. [15]. Being Extended to finite temperature, the effective two-body 
interaction V^^lT) in r-space is given by, 

{f,'f,'\Vf{T)\m) = iTr^/(fc„,p,/?,T) J dr3df3>;(r3')(l-r/(ri3,T))(l-r7(r^3,r)) 

X W3{f^'f2'f^'\fif2f3)(pn{r3){l - r]{ri3,T)){l - ?7(r23,r)) (7) 

where the trace is taken with respect to the spin and isospin of the third nucleon. The func- 
tion r){r,T) is the correlation function ( or called defect function ) [25]. The defect function 
is directly related to the solution of the finite-temperature Bethe-Goldstone equation, i.e. 
eq.(l) [23,25] and depends on temperature. The transformation of the TBF to the above 
effective interaction entails a self-consistent coupling between the TBF and the Brueckner 
procedure of solving the Bethe-Goldstone equation. One first calculates the correlation func- 
tion with only the two-body force and then builds up the effective interaction V^'^'^(r) which 
in turn is added to the bare two-body force, and again evaluates the correlation function 
and so on up to convergence is reached. 

It is worth to noticing that the TBF itself is the same as that adopted in the 
zero-temperature case [17] and independent of temperature. As compared to the zero- 
temperature case, the effective two-body force V^^{T) constructed from the TBF depends 
on temperature due to the medium effects and consequently its effects are expected to be 
more pronounced at finite temperature. It is obviously from Eq.(7) that the temperature 
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dependence of Vf{T) comes from the Fermi distribution f{k, p, {3, T) and the defect func- 
tion ri{r^ T) . This temperature dependence is somewhat similar to that of the Skyrme force 
where a three-body part is equivalent to a two-body density-dependent effective interaction 
which depends implicitly on temperature via density [11]. 

Wc solve Eqs.(l), (4), (5), (6), (7) self-consistently to get the G-matrix for any given 
density, temperature and isospin asymmetry. In general, five iterations is needed to reach a 
satisfactory convergence. The total energy E{p,P,T) — Ekin{p, (3,T) + Epot{p, P,T) can be 
obtained from the G-matrix. The total kinetic energy is 

EUp, /3, T) = ^ ^ U{k, p, (5, T)^ (8) 

T k 

and the total potential energy 

£;,„i(p,/9,T) = l^^/.(fc,p,/3,T)/,,(fc',p,/9,T)(A;A;' | G.,.,(p, A T, e(fc) + e(^')) I kk')^ (9) 

t,t' k,k' 

The total entropy S can be evaluated in the approximation of a non-interacting Fermi gas 
of quasi-particle in the mean field Ur{k, p, [3, T) [14]. The free energy and the pressure can 
be obtained according to the standard thermodynamic relation, i.e, 

F^E-TS (10) 

Based on the above extended FTBHF approach including the TBF, we can investigate the 
thermodynamic properties of hot asymmetric nuclear matter at various values of density, 
temperature and asymmetry. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The EOS of symmetric nuclear matter at finite temperature is reported in Fig. 2, 
where the isotherms of pressure versus density is given for six values of temperature 
T = 0, 8, 10, 12, 14, 16MeV from the bottom to the top. In the figure the solid and dashed 
curves indicate the results obtained by using the AVig two-body force plus the TBF and 
the pure AVis two-body force, respectively. As expected, the TBF gives a repulsive con- 
tribution to the EOS of nuclear matter. This contribution becomes stronger as increasing 
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density and makes the EOS at high density much stiffer. At zero temperature, the additional 
repulsion from the TBF leads to a decisive improvement of the calculated saturation density 
of symmetric nuclear matter towards the empirical value ( comparing the lowest solid and 
dashed curves ). It is clearly seen from the figure that in both cases with and without 
the TBF, the predicted EOS exhibits a typical Van der Waals structure in agreement with 
the results from the Skyrmc-Hartrcc-Fock calculations [9] and from the FTBHF approach 
based on two-body NN interactions [26], which indicates that the infinite nuclear matter 
can undergo a liquid-gas phase transition. Without the TBF, the critical temperature Tq 
of the hquid-gas transition is approximately 16MeV. Inclusion of the TBF reduces the crit- 
ical temperature to about 13MeV which is smaller than the values Tc — 15 — 20MeV from 
the Skryme-Hartree-Fock calculations [9]. The TBF reduction of the critical temperature 
is readily understood since the TBF effect becomes more pronounced as increasing density 
and temperature. In finite nuclei, the inclusion of the finite size effects and the Coulomb 
force may lead to a considerable reduction of Tc as discussed in Rcf. [11]. The present value 
Tc — 13MeV is comparable with the value Tc — 14MeV predicted by the relativistic mean 
field theory [5] but it is larger than the ones obtained from the Dirac-Brueckner ( DB ) 
approach, 7b ^ 8 ~ 9MeV in Ref. [12] and Tc ~ lOMeV in Ref. [13]. One possible reason 
for this discrepancy is the difference between the relativistic effect in the DB approach [27] 
and the effect of the present microscopic TBF as discussed in Ref. [17]. 

In Fig. 3 is reported the proton or neutron s.p. potential in symmetric nuclear matter 
at a fixed density p = 0.16fm^^ for three different values of temperature T = 0, 10 and 
20MeV. As increasing temperature, the s.p. potential becomes more repulsive, resulting in 
an enhancement of the potential part of the energy per nuclcon in nuclear matter at finite 
temperature. It is also seen that the curvature around the Fermi momentum becomes more 
smooth as increasing temperature. This is in agreement with the prediction obtained by 
using only pure two-body forces [26] and is attributed to the thermal excitations around the 
Fermi surface at finite temperature. 

Due to the isospin effect, the proton and neutron s.p. potentials are different in asym- 
metric nuclear matter. In order to discuss the TBF effects on the isospin dependence of 
the proton and neutron s.p. potentials, in Fig. 4 is depicted the proton and neutron s.p. 
potentials at momentum A; = as a function of asymmetry parameter for p = 0.16fm~^ 
and T — lOMeV. In the figure, the solid and dashed curves indicate the results obtained by 
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adopting the AVis plus the TBF and the pure AVis two-body force, respectively. It is clear 
from the figure that the inclusion of the TBF in the calculations leads to an additional repul- 
sion for both the proton and neutron s.p. potentials in the whole isospin asymmetry range 
< /? < 1. As similar to the zero-temperature case, the proton s.p. potential becomes more 
attractive while the neutron one becomes more repulsive going from symmetric (/3 = 0) to 
neutron (/9 = 1) matter. Such a different isospin dependence of the proton and neutron s.p. 
potentials stems mainly from the attractive contribution of the isospin singlet SD tensor 
channel which becomes stronger for protons and weaker for neutrons as increasing neutron 
excess [22]. As compared to the two-body force predictions, the TBF contribution shift 
the proton and neutron s.p. potentials to higher values. In both cases with and without 
the TBF contribution, the dependence of the proton and neutron s.p. potentials on the 
asymmetry parameter /3 is almost linear as the same as in the zero-temperature case, which 
supports microscopically the Lane assumption [28] and extends its validity to the case of 
finite temperature. 

In heavy-ion collisions using neutron-rich nuclei, hot nuclear matter of high isospin asym- 
metry could be formed and the reaction dynamics is expected to be very sensitive to the 
density and temperature dependence of the symmetry energy. In the prompt explosion 
model of a type-II supernova the /3-capture process drives the star to the last stage of the 
collapse where the proton fraction is about 0.33 and the temperature of a few tens of MeV 
may be reached. Therefore, the understanding of the cooling mechanism of a protoneu- 
tron star requires the information of the symmetry energy in hot nuclear matter with high 
accuracy [29,30]. 

In Fig. 5 is depicted the shift of the energy per nuclcon between asymmetric nuclear 
matter and symmetric nuclear matter, Ea{p, P,T) — Ea{p, P — 0,T), versus (5^ at T — 
20MeV for two different densities p — 0.16 and 0.32fm~^. It is shown that the energy per 
nucleon of asymmetric nuclear matter at finite temperature fulfills a quadratic dependence 
on asymmetry parameter /3 in the whole asymmetry range ( < /3 < 1 ) as the same as 
in the case for zero temperature [22,31,32]. The above /3^-law comes originally from the 
empirical mass formula for small isospin asymmetry. The present calculations show that the 
validity of the empirical /9^-law can be extended to the highest isospin asymmetry and to 
the case of finite temperature. 

Besides the /3^-dependent term, an additional binding term called Wigner energy has 
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been found for the N — Z nuclei [33] . The Wigner energy is hnear in (3 and manifests itself 
as a spike in the experimental isobaric mass parabola. The Wigner energy has been explored 
and discussed extensively in the literatures [34]. It is generally expected that the Wigner 
energy originates from the neutron-proton {np) pairing correlations and is the strongest 
prefigure of the np paring in the N = Z nuclei [35]. In the present calculations, we do not 
consider the influence of the np pairing correlations. Since the np pairing goes to vanishing 
very rapidly as soon as the isospin asymmetry and the temperature deviates from zero [36], 
we expect that the effect of the np pairing correlation shall be appreciable only for the cold 
symmetric nuclear matter, which is found to slightly lower the energy per nucleon of the 
cold symmetric nuclear matter by a self-consistent treatment of the BCS gap equation and 
the BHF equations [37]. 

The symmetry energy is defined as 

Due to the simple quadratic dependence of the energy per nucleon on asymmetry parameter, 
the symmetry energy can be equivalently calculated as the difference between the energy 
per nucleon of pure neutron matter and symmetric nuclear matter, i.e., 

EsymiP,T) = EA{p,f3^1,T)-EA{p,P^0,T). (13) 

To see the infiuence of temperature on the symmetry energy, in the left panel of Fig. 6 we 
show the density dependence of the symmetry energy for three different temperatures T — 
0, 10 and 20MeV. It is seen that the symmetry energy decreases as increasing temperature 
and at finite temperature it remains a monotonic increasing function of density. As compared 
to the result by using only a two-body force [38], the inclusion of the TBF contribution 
enlarges the sensitivity of the symmetry energy with respect to the variation of temperature, 
especially at high densities. For a given temperature, the TBF contribution leads to a strong 
enhancement of the symmetry energy at high densities as shown in the right panel of Fig. 6 
where the results for T — 20MeV obtained by using the AVis plus the TBF ( sohd curve ) 
and without the TBF ( dashed curve ) are plotted. It is also seen that the symmetry energy 
predicted by including the TBF contribution rises much more steeply with density than the 
corresponding two-body force prediction. 

Before summary, we discuss briefly the temperature dependence of the effective mass and 



d^EA{p,P,T) 



(12) 
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the TBF effect on the effective mass. The effective mass m* for neutron (r = n) or proton 
T = p is defined as [14] 



where er{k) is the neutron or proton s.p. energy. The effective mass describes the nonlocal 
part of the nuclear mean field which makes the local part less attractive for a nucleon 
travelling with momentum A; > 0. In Fig. 7 is shown the momentum dependence of m* 
for symmetric nuclear matter at a fixed density p — 0.16fm~^ and three different values of 
temperature T — 0,10 and 20MeV. The results are comparable with the ones reported in 
Refs. [39,40]. As similar to the zero-temperature case, the momentum dependence of m*{k) 
is characterized by a wide bump around the Fermi momentum due to the high probability for 
particle- hole excitations nearby Fermi surface [23]. However as the temperature increases, 
the peak of m* becomes flatter and the peak value of m* becomes lower. These results 
are in agreement with the previous calculations by adopting only two-body forces [14,25,26] 
and related directly with the temperature effect on the s.p. potential around the Fermi 
momentum. 

To see the isospin effects, in Fig. 8 we display the neutron and proton effective masses 
at their respective Fermi momenta kp and kp for T — lOMeV as a function of asymmetry 
parameter f3. In the figure, the sohd curves are calculated by using the AVis plus the TBF, 
while the dashed curves by adopting only the AVis two-body force. In both cases with and 
without including the TBF contribution, the neutron effective mass m* increases and the 
proton one m* decreases as increasing /3. The TBF effect is to shift m* and m* to slightly 
higher values. One can also see from the figure that the linear scissor-shaped behavior 
predicted for T = [22] remains unchanged in the case for finite temperature. 



In summary, we have introduced the TBF contribution into the FTBHF approach and 
investigated the EOS of asymmetric nuclear matter at finite temperature. By using the 
extended model, we find that the TBF contribution to the EOS is repulsive. The TBF 
repulsion increases rapidly as increasing density and consequently the high-density EOS 
including the TBF contribution becomes much stiffer as compared to the one obtained by 
adopting only the two-body force. Within the extended model, the calculated EOS of hot 




(14) 



IV. SUMMARY AND CONCLUSION 
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symmetric nuclear matter exhibits a typical Van der Waals structure, implying the presence 
of a liquid-gas phase transition. The obtained critical temperature of the liquid-gas phase 
transition is roughly 16MaV in the case without including the TBF contribution. The TBF 
leads to a 3MeV reduction of the critical temperature to about 13MeV which is between 
Tc = 15 — 20MeV from the non-relativistic approaches and Tc — lOMev from the Dirac- 
Brucckner method. 

Our calculations show that the simple /3^-law fulfilled by the energy per nucleon of 
asymmetric nuclear matter at T = can be extended to the case of finite temperature. This 
is of interest, since it means that the energy per nucleon of hot asymmetric nuclear matter 
can be extracted from the two limiting cases of symmetric nuclear matter and pure neutron 
matter. The symmetry energy has been calculated. For a given temperature, the symmetry 
energy is a monotonic increasing function of density in both cases with and without the 
TBF. At a fixed density, the symmetry energy decreases as the temperature increases. The 
TBF gives a strong enhancement of the symmetry energy at high density and makes the 
symmetry energy rise much more rapidly as compared to the results without including the 
TBF contribution. 

The temperature dependence and the isospin dependence of the proton and neutron s.p. 
potentials and effective masses have also been calculated and discussed. As expected, the 
s.p. potential becomes more repulsive and the peak of the effective mass around the Fermi 
momentum becomes flatter as increasing temperature. It is shown that the neutron effective 
mass increases and the proton one decreases as increasing /?. Both the TBF contribution 
and the finite temperature effect do not change this isospin behavior of the neutron and 
proton effective masses. 

Our results show that the TBF contribution affects considerably the stiffness of the EOS 
of hot nuclear matter and the symmetry energy at high density, and as a consequence, it 
may have important implication for the investigations of the dynamics evolution of supernova 
explosions and the structure of protoneutron stars. 
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FIG. 1. 
Ref.[15] 



FIGURES 

Diagrams of the microscopic TBF adopted for the present calculations, taken from 



FIG. 2. Pressure as function of density for symmetric matter at six values of temperature 
r = 0, 8, 10, 12, 14, 16MeV from the bottom to the top. The solid and dashed curves are the results 

by using the ^Vig plus the TBF and the pure AVis two-body force, respectively 

FIG. 3. Momentum dependence of the single-particle potential in symmetric nuclear matter at a 
fixed density p = 0.16/m~^ and three different temperatures T = 0, 10, 20MeV. In the calculations, 
the TBF is included. 

FIG. 4. Proton and neutron single-particle potentials at = versus asymmetry parameter 
/? for /) = O.iefm"''^ and T = lOMeV. The solid curves arc the results by including the TBF 
contribution and the dashed ones without the TBF contribution. 

FIG. 5. Energy per nucleon Ea{p, P, T) — Ea{p, P = 0,T) versus 0^ in the range < /3 < 1 for 
two different densities. The results are obtained by including the TBF. 

FIG. 6. Left panel: Symmetry energy as a function of density for three values of temperature. 
The results are obtained by including the TBF. Right panel: Density dependence of the symmetry 
energy with the TBF included (solid line) and not included (dashed line) at a fixed temperature 
T = lOMeV. 

FIG. 7. Momentum dependence of the nucleon effective mass in symmetric nuclear matter at 
three different values of temperature. 

FIG. 8. Proton and neutron effective masses as a function of asymmetry parameter obtained 
by adopting the AVis plus the TBF ( solid curves ) and by using only the AVis two-body force ( 
dashed curves ). 
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